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Abstract. We establish several new stylised facts concerning the intra-day seasonalities of stock 
dynamics. Beyond the well known U-shaped pattern of the volatility, we find that the average 
*^ correlation between stocks increases throughout the day, leading to a smaller relative dispersion 

I between stocks. Somewhat paradoxically, the kurtosis (a measure of volatility surprises) reaches a 

minimum at the open of the market, when the volatility is at its peak. We confirm that the dispersion 
^■vq kurtosis is a markedly decreasing function of the index return. This means that during large market 

swings, the idiosyncratic component of the stock dynamics becomes sub-dominant. In a nutshell, 
early hours of trading are dominated by idiosyncratic or sector specific effects with little surprises, 
whereas the influence of the market factor increases throughout the day, and surprises become more 
frequent. 

j 1 1. Introduction 

m Financial markets operate in sync with human activities. It is therefore no surprise that financial 

time series reveal a number of seasonalities related to human rhythms: markets open in the morning 
and close in the evening, remain closed during week-ends and during vacations; wages are paid and 
portfolios are re-balanced on a monthly basis, earnings are announced on a quarterly basis (in the US), 
i ^7 etc. These periodicities leave a statistical trace on the time series of returns of many assets. Among 

the best known periodicities is the so-called U effect [TJ that describes the intra-day pattern of 
volatility of individual US stocks: the average volatility is observed to be high after the market opens, 
then decreases as to reach a minimum around lunch time and increases again steadily until market 
close0 

[--^ In this short note we want to report on additional intra-day patterns concerning both individual and 

collective stock dynamics. First, we study the intra-day pattern of other moments of the individual 
stock dynamics, beyond the well known U-shaped volatility. Second, we characterise the cross-sectional 
distribution of returns and its typical evolution during the day. Finally, we study the correlation 
matrix between stock returns and find that the leading modes also have a very well defined intra-day 
pattern. Our study here is entirely empirical, but our results certainly beg for a detailed theoretical 
^ interpretation in terms of agent behaviour: strategies, information processing, risk aversion, etc. We 

J^~j provide some hints in that direction in the conclusion. 

rS 

2. Data, notations and definitions 

We have considered a set of N = 126 stocks of the New York Stock Exchange (which are among 
the 250 largest market capitalisations) that has been continuously traded during the period between 
01/01/2000 and 12/31/2009 to form a statistical ensemble of 5 minutes stock returns. The total 
number of 5 minute bins is 186, 498, whereas the total number of days is T — 2, 391 (K = 78 bins per 
day). Stocks will be labelled by a = 1, . . . , N, days by t = 1, . . . , T and bins by k = 1, . . . , K . The 
return of stock a in bin k of day t will be denoted as rj a (k] t). 

Different types of averages will be needed. Time averages for a given stock and a given bin are 
expressed with angled brackets: (...), whereas averages over the ensemble of stocks for a given bin 
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Date: September 27, 2010. 

^This pattern is a little different in Europe or in the UK, with a second volatility spike at 2:30 pm GMT when the 
US market opens. 
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in a given day appear with square brackets: [...]. For an arbitrary function F(.) of these returns we 
therefore write: 

T NT 

(F)(k;a):=^J2F(r, a (k;t)); [F](k;t) := 1 £ F( Va (k;t)); [(F)](k) = ^J2[F](k;t). (1) 

t=l a=l t=l 

The first set of observables concerns single stock properties. We characterise the distribution of stock 
a in bin k by its four first moments: mean /j, a (k), standard deviation (volatility) o~ a (k), skewness £ Q (fc) 
and kurtosis n a {k). We will in fact use low moments, less noisy estimates of the last two quantities. 
We will define m a {k) as the median of all returns of stock a in bin k, and define: 

Ha(k) = (ri a (k;t)) (2a) 

a 2 a (k) = {Va(k;t) 2 ) - (2b) 

Ca(k) = — 7TT (Ma(fc) ~ m a (fc)) (2c) 

KqW = 24 v V 2 >J +<«(*) ' (2d) 

Within a cumulant expansion, the last two lines coincide with the usual definition of skewness and 
kurtosis, but no moments larger than two are needed to estimate them. Note that the correction term 
Ca(k) 2 to the kurtosis turns out to be negligible, and we have neglected it in the following. We will 
be interested below in the average over all stocks of the above quantities, as a way to characterise the 
typical intra-day evolution of the distribution of single stock returns. 

One can also consider cross sectional distributions, i.e. the dispersion of the returns of the N stocks 
for a given bin k in a given day t, i.e. one distribution every five minutes. One can again characterise 
these distributions in terms of the first four moments. The median of all N returns for a given k; t is 
now uid(k;t) (d for "dispersion"), and we define: 

Hd{k:t) = [r) a (k;t)] (3a) 
a 2 d (k;t) = [r ]a {k;t) 2 ]- l i 2 d {k;t) ( 3b ) 

Cd(M) = — tt-ty (Vd(k;t) - m d (k;t)) (3c) 

Note that fid(k;t) can be seen as the return of an index, equiweighted on all stocks. We will be 
interested below in the average over all days of the above quantities, as a way to characterise the 
typical intra-day evolution of the dispersion between stock returns. 

Although the dispersion already captures part of the "Co-movements" of stocks, a more direct 
characterisation is through the standard correlation of returns. In order to measure the correlation 
matrix of the returns, we first normalise each return by the dispersion of the corresponding bin. 
This factors in any "trivial" intra-day seasonality, and also accounts for the fact that the volatility 
fluctuates quite a bit during the 10 year time interval that we consider. Therefore, we introduce: 
rf a (k;t) — rj a (k]t)/ad(k;t) and study the correlation matrix defined for a given bin k: 

c*Ak) ■= zjj^-jk) (fr«(M)*fc(M))c, (4) 

where the subscript c means "connected part" (i.e. averages have been subtracted) and := 
(?7 2 (fc;t)) c . Of special interest are the largest eigenvalues and eigenvectors of C a ,p(k), which chara- 
terize the correlation structure of stock returns. This analysis has been performed in several papers 
(0 [3 I17[ I14j ) using daily or high frequency returns, and it is well known that the structure of large 
eigenvectors reflects the existence of economic sectors of activity. The largest eigenvalue Ai , in partic- 
ular, corresponds to the market mode, and is associated to an eigenvector with all entries positive and 
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Figure 1. We show the average volatility of stocks cr(fc), the average cross sectional 
dispersion Gd{k) and the average absolute value of the index return {\jjL,i(k,t)\) (mul- 
tiplied by 4 for clarity) as a function of k, with the corresponding statistical error 
bars. All display the well known U pattern. We also compare a(k) with a power-law 
decay with j3 sa 0.3 (dashed line), which is a good fit for the first half of the day. 
Inset: ratio a(k)/o~d(k) as a function of fc, showing that dispersion effects diminish 
throughout the day. 

close to 1/^/N. In fact, Xi/N can be seen as a measure of the average correlation between stocks. We 
will be interested below in the k dependence of the largest eigenvalues and their associated eigenvectors, 
a study that, to the best of our knowledge, has not been reported in the literature before. 

3. Single stock intra-day seasonalities 

3.1. Odd moments. Odd moments tend to be small and noisy, so it is difficult to draw definite 
conclusions. The average return is on average over the whole period positive, but noisy and does 
not show any intra-day pattern. The average skewness of five minutes returns is also noisy and is 
compatible with zero, again without any identifiable intra-day pattern at all. This is at variance with 
the skewness of returns on a longer time interval, which is negative. The build up of negative skewness 
with time scale is a consequence of the leverage effect, i.e. negative returns tend to be followed by 
larger volatilities (see e.g. [H]). 

3.2. Even moments. The average volatility, on the other hand, reveals a very clean U-shaped pattern 
that has been reported many times in the literature ([TJI2])- We show in Fig. [I](7(fc) = [a a (k)]. Note 
that the overnight volatility (=s 1.15% is much larger than the typical five minute volatility, and is not 
shown in the graph. Interestingly, the average volatility is found to decay in the first two hours of 
trading as a power-law with j3 sa 0.3. This relaxation is reminiscent of the power-law decay of 
the volatility after large price swings [TT1 [T51 [TBI [TO] . The overnight return is indeed usually quite 
large, and can be seen as a strong perturbation. The power-law relaxation suggests that some critical 
mechanism is involved in the way volatility reverts back to 'normal' after market jumps. 

Turning now to the kurtosis n(k) — [k q (/c)], we find, perhaps surprisingly, that there is a clear intra- 
day growth of the kurtosis from k sa 3.5 at the beginning of the day to k sa 5 around 1 p.m., and stays 
approximately constant (but noisy) until the end of the day - see Fig. [5] left. The overnight kurtosis 
remains around 5. This finding is counter-intuitive because one would naively associate the large 
volatility in the morning with huge swings, symptomatic of the market uncertainty at the open. But 
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FIGURE 2. Left: average kurtosis of individual stocks n(k) as a function of k. Right: 
average cross sectional kurtosis Kd{k) as a function of k. In both cases, we show the 
l-o" dispersion around the mean (i.e. not the error bar). The red triangles correspond 
to the overnight values. 



this is not the case: the maximum of the intra-day volatility corresponds to a minimum in kurtosis. 
We will report similar counter-intuititive results below. Possible mechanisms are discussed in the 
conclusion. 

4. Cross-sectional intra-day seasonalities 

4.1. Odd moments. Noting that the average over stocks of fj, a (k) is identical to the average over 
time of /id(fc; t), the discussion of the first moment of the cross-sectional distribution is redundant. The 
average of |/Xd(fc;t)| is a proxy for the index volatility, and is displayed in Fig. [I]: it shows a U-shaped 
pattern similar to that of cr(k), with however a stronger end-of-day surge. This is due to the correlation 
pattern discussed in section 5 below: the average correlation between stock indeed increases as the day 
proceeds, leading to an increased index volatility. 

As far as the average skewness Cd(k) = (Q(fc;i)) is concerned, we again find a very noisy quantity 
with no particular intra-day pattern. The only notable feature is that this time, the skewness is 
significantly positive, albeit small: the average over k of (d(k) is found to be ~ 0.025. 

4.2. Even moments. As above, the even moments show clear patterns. The average dispersion 
Cd(fc) = (&d(k;t)) exhibits a U-shaped pattern very similar to that of cr(fc) — see Fig. [I] In fact, 
the ratio a{k)/cjd(k) is plotted in the inset of Fig. [l]as a function of k and increases from rj 1.15 at 
the open to 1.45 at the close. In relative terms, the dispersion is thus stronger in the morning, and 
decreases as the day proceeds. 

The dispersion kurtosis Kd(k) = (Kd(k;t)}, on the other hand, has an invertedV shape, and reaches 
a minimum at the open and at the close of the market, i.e. when the dispersion and the volatility are 
locally maximum. So even when the dispersion of returns is at its peak, with stocks all over the place 
(so to say), the cross-sectional distribution of returns is on average closer to a Gaussian! Note however 
that the variation of the kurtosis is not large, from Kd = 2 to Kd = 2.4. The overnight dispersion 
kurtosis, on the other hand, is much stronger: Kd ~ 3.3. 
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Figure 3. Cross sectional dispersion ad as a function of the index return (equi- 
weighted on all stocks) fid- We added error bars, and linear branches that fit the 
small \fj,d\ slopes, that emphasise the sub-linear behaviour of ad- 
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FIGURE 4. Cross sectional skewness (d as a function of the index return (equiweighted 
on all stocks) fid- We added error bars (dotted lines), that are actually difficult to see 
near the origin. 

4.3. Conditioning on the index return. As noted above, the quantity fid{k\t) is the return of 
an equiweighted index. It is interesting to condition the value of the moments of the cross-sectional 
dispersion on this quantity. Such a study was performed on daily returns in [13] and more recently by 
L. Borland 0]. In agreement with the results of [21 E], we find that the average dispersion ad is an 
increasing function of the amplitude of the index return, see Fig. [3} As noted in [9], this observation 
shows that the volatility of the stock residuals in a one-factor model must depend on the volatility 
of the market mode. Fig. [3] furthermore suggests that this dependence is sub-linear (see [S] for some 
elaborations on this observation.) 
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Figure 5. Left: Cross sectional kurtosis function of the index return (equi- 

weighted on all stocks) fid- Note that the dependence is nearly the same for positive 
and negative market returns. Right: Cross sectional kurtosis Kid clS el function of the 
cross sectional dispersion er^. We added error bars (dashed lines) on both figures. 

As first established in [12] on daily data, we find that the skewness Cd is an odd function of fid, as 
shown in Fig. |4j Note that the skewness increases very abruptly for small fid and saturates for larger 
values of the index return. Pictorially, a positive index return can be thought of as resulting from a few 
"winners" running ahead of the pack, contributing both to the mean fid and to the skewness. The slope 
of Cd((J>d) around the origin, together with the fact that the index has made on average positive daily 
gains in the period 2000 - 2009, are enough to explain the average value of the dispersion skewness 
Cd(k) ~ 0.025 reported above. 

Finally, the dispersion kurtosis Kd shows again a non-intuitive decreasing behaviour as a function of 
\(id\, see Fig. [5} The average kurtosis conditioned to a value of \fid\ decreases from w 2.8 for small index 
returns to ~ 1.8 for index returns larger than 2% in absolute value. This was first noticed en passant 
in [9 on daily data and recently emphasised by Borland [4]. Here, we confirm on five minute returns 
this strange stylised fact: the cross-sectional distribution of returns appears to be more Gaussian when 
its mean is off-centred. 

However, if we now condition Kd on the dispersion ad (which, as we found above, is positively 
correlated with \fid\), we find (see Fig. [5]) the opposite behaviour, i.e. the larger the dispersion, the 
larger the kurtosis k^! We will offer a discussion of these confusing effects in the discussion section 
below. 

5. Intra-day seasonalities in the inter-stock correlations 

Let us now turn to the properties of the eigenvalues and eigenvectors of the N x N correlation 
matrix C at p(k) defined by Eq. Q above. 

5.1. The top eigenvalue. The largest eigenvalue Ai of the correlation matrix of stock returns is well 
known to be associated with the "market mode" , i.e. all stocks moving more or less in sync. As 
recalled above, the quantity Xi/N can be used to define the average correlation between stocks. 

We show in Fig. [6] (left side) the magnitude of X±/N as a function of k. Interestingly, the average 
correlation clearly increases as time elapses, from a rather small value ss 0.12 when the market opens 
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FIGURE 6. Left: Dependence of the top eigenvalue of the correlation matrix C(k), 
Xi(k)/N, as a function of time of day. Right: Evolution of the scalar product between 
top eigenvector vi(k) and the uniform vector e. The red triangles correspond to the 
overnight values. 

to w 0.3 near market close. This is in agreement with the fact that the dispersion ad(k) is, in relative 
terms, smaller at the end of the day (see Fig. [I] inset). The value of Xi/N for the correlation of 
overnight returns is also around 0.3, in continuity with the value at the end of the trading day. 

In agreement with the idea that the stock dynamics become more and more uniform as the day 
proceeds, we find a substantial increase of the scalar product of the largest eigenvector v±(k) with the 
uniform normalised vector e = (l/*/N, l/*/N, . . . , l/\/jV) — see Fig. fright. This scalar product 
is always close to unity, confirming the market mode interpretation of the top eigenvalue, but starts 
the day around 0.97 and ends the day at 0.995, before dropping again in the last bins of the day and 
during the overnight, when it is equal to 0.985 (i.e. larger than the open value). 

5.2. Smaller eigenvalues. The evolution of the next six eigenvalues Aj(fc), i = 2, . . . , 7 is shown in 
Fig. [7J We see that the amplitude of this risk factors now decreases with time, before shooting back 
up during the overnight (see the last point of the graphs). Although by construction the trace of 
the correlation matrix, and therefore the sum of all N eigenvalues is constant (and equal to N), this 
decrease is not a trivial consequence of the increase of Ai, since the sum of the first five eigenvalues 
is ~ 50, still small compared to Tr(C) = N = 126. What we see here is that as the day proceeds, 
more and more risk is carried by the market factor, while the amplitude of sectorial moves shrivels in 
relative terms (but remember that the correlation matrix is defined after normalising the returns by 
the local volatility, which increases in the last hours of the day). 

It is more difficult to visualise the evolution of the corresponding eigenvectors, since there is no 
natural vectors to compare them with. Furthermore, eigenvalues can "collide" and cross, resulting in 
an interchange between two consecutive eigenvectors. We have therefore chosen to take as a reference 
the eigenvectors Vi(l) in the opening bin k = 1, corresponding to the largest values of Aj(fc), i = 2, . . . , 7. 
We then form the 6x6 matrix of scalar products Wij(k) = Vi(l)-Vj(k). The singular values si(k) of this 
matrix (equal to the square-root of the eigenvalues of W T W) give a measure of the overlap between 
the eigenspace spanned by the Vi(l) and that spanned by the vAk). If the vj(k) are a permutation 
of the Uj(l), all the Sj's are equal to unity, indicating maximum overlap. In particular, se(l) = 1 
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FIGURE 7. Left: Smaller eigenvalues Xi(k), i — 2, ... ,7) as a function of k. Right: 
The 6 singular values sg(k) (£ — 2, ... ,7) of the matrix Wij(k). In the absence of 
any true persistence, one would expect all singular value to lie in the interval [0, 0.12], 
much below the smallest singular value S7 ~ 0.8. The triangles correspond to the 
overnight values. 

trivially. The evolution of the Sf(k), £ — 2, ... ,7 is shown in Fig. [7] right. Using the results of [5], we 
conclude that all si(k) are meaningful, since in the absence of any true correlations between the Vi(l) 
and the Vj(k), one would expect all singular values to lie in the interval [0,0.12]. Therefore, although 
the structure of correlations clearly evolves between the opening hours and the closing hours, there is 
as expected a strong overlap between the principal components throughout the day. 

6. Discussion & Conclusion 

Let us present a synthetic account of the above empirical results, for which we only propose an 
interpretation stub. We have seen that during the opening hours of the market, the volatility and 
the dispersion of returns are high, whereas kurtosis effects are relatively low. These two quantities 
are different measures of the heterogeneity of stock returns, and quite paradoxically they are found 
to behave in opposite ways. But while the volatility and dispersion are dimensional measures of 
heterogeneity (measuring the spread of returns in %), the kurtosis is a relative, a-dimensional measure 
of surprise. What our results mean in intuitive terms is that although the typical amplitudes of stock 
returns are high in the morning, outliers are relatively rare, both over time and over stocks. In a sense, 
agitation is the norm during these early hours of trading, stocks move in different directions in such 
a way that the average correlation is weaker than average, and the top eigenvector of the correlation 
matrix is farther away from the uniform mode e = (1/vJV, 1/y/N, . . . , 1/yN). But anomalously large 
jumps rarely take place in the morning — as expected, these jumps are more likely overnight (and are 
to be related to arrival of corporate specific or market-wide information), where kurtosis effects are 
strongest, both for single stock and cross-sectional returns. As the day proceeds, correlations increase 
and dispersion decreases, but unexpectedly large jumps become more probable, thereby increasing the 
kurtosis. 

The second somewhat paradoxical effect is the dependence of the kurtosis on the index return, 
which was recently interpreted by L. Borland as a signature of collective behaviour during crises 0]. 
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Again, days when the market as a whole moves a lot are also large dispersion days where all stocks 
move a lot in different directions, but with little outliers, i.e. one or a handful of stocks that would 
jump up or down. In this sense, these days are more homogeneous. Should one deduce from this that 
there is a stronger "synchronisation" , or collective dynamics, during these periods, as suggested by 
Borland? While it is true that the average correlation between stocks depends on the index return, 
this dependence is in fact signed: correlations are stronger for negative index returns and weaker 
for positive returns, see [21 115] . This is in contrast with the kurtosis effect discussed here, which is 
surprisingly symmetrical (see Fig. [5|. A quantitative model for this behaviour is missing at this stage. 
Qualitatively, however, we believe that the mechanism is the following [5]: when the index return 
is large, the dominant source of dispersion becomes the market exposure (the '/3's') of the different 
stocks, rather than the idiosyncratic residuals. Since the distribution of the /3's is roughly Gaussian, 
kurtosis effects do indeed decrease for large index returns. This interpretation however requires that 
the volatility of the residuals increases sub-linearly with the index volatility, as indeed suggested by 
the data shown in Fig. [3] The fact that during large swings of the index, the market exposure of stocks 
becomes the dominant factor is probably a result of index/futures arbitrage. 

Finally, although large index return days are large dispersion days, the converse is not true. A 
typical large dispersion day is in fact a day when one or a handful of stocks gyrate wildly, contributing 
both to the dispersion and to the kurtosis, and explaining the positive correlation between er^ and 
Kd- If this interpretation is correct, this positive correlation should diminish when one uses the mean- 
absolute deviation and not the variance to compute the dispersion, since the former is less sensitive to 
outliers. We have checked that this is indeed the case. 

To summarise, we have established several new stylised facts concerning the intra-day seasonalities 
of stock dynamics. Beyond the well known U-shapcd pattern of the volatility, we have found that 
the average correlation between stocks increases throughout the day, leading to a smaller dispersion 
between stocks (in relative terms). However, the kurtosis, which is a measure of volatility surprises, is 
in fact minimum at the open of the market, when the volatility is at its peak. We have also confirmed 
that the dispersion kurtosis is a symmetric, markedly decreasing function of the index return. This 
means that during large market swings, the idiosyncratic component of the stock dynamics becomes 
sub-dominant, an effect that we have confirmed directly. Finally, while the market mode component 
of the dynamics becomes stronger as the day proceeds, the sectorial components recede. In a nutshell, 
early hours of trading are dominated by idiosyncratic or sector specific effects with little surprises, 
whereas the influence of macro, market factor increases throughout the day, and surprises become 
more frequent. A detailed quantitative interpretation of our results, for example of the power-law 
decay of the volatility in the morning, is at this stage lacking. We believe that, when available, such 
an interpretation will shed light on the relative importance of behavioural and informational effects on 
price formation and volatility. 
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